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ABSTRACT 

A relativistic colhsionless shock propagating into an unmagnetized medium leaves 
^^ ■ behind a strong large-scale magnetic field. This seems to follow from two assumptions: 

^Sj . (i) GRB afterglows are explained by synchrotron emission of a relativistic shock, (ii) 



^ 



magnetic field can't exist on microscopic scales only, it would decay by phase space 
mixing. Assumption (i) is generally accepted because of an apparent success of the shock 
l/~) [ synchrotron phenomenological model of GRB afterglow. Assumption (ii) is confirmed 

in this work by a low-dimensional numerical simulation. One may hypothesize that 
^ ■ relativistic shock velocities are not essential for the magnetic field generation, and that 

all colhsionless shocks propagating into an unmagnetized medium generate strong large- 
scale magnetic fields. If this hypothesis is true, the first cosmical magnetic fields could 
have been generated in shocks of the first virialized objects. 
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1. GRB Afterglow Phenomenology 

It appears that GRBs and GRB afterglows can be adequately described by synchrotron emis- 
sion of relativistic shocks (eg. Waxman 1997, Sari et al 1997, Guetta et al 2000). These phenomeno- 
;h ' logical models describe the shock by a set of adjustable parameters. One of these parameters, 

^B = -B^/(87re), is the magnetic energy density divided by the proton energy density e downstream 
of the shock transition. It is thought that ^b is close to unity, say ^^ ~ 0.1 . Much smaller values, 
(,B ~ 0.001, were given by Panaitescu & Kumar (2000), but as we are about to show even 0.001 is 
a surprisingly large number. 

For the sake of simplicity, consider a mildly relativistic stage of a GRB afterglow, when Lorentz 
factor is ~few. Assume that unshocked plasma is at a typical ISM temperature, ~ leV. The 
unshocked magnetic energy per particle is ^ leV, because the field has to be confined by the 
plasma pressure. If hydrodynamics is applicable, then the blast wave is stable (Gruzinov 2000), 
and there is no magnetic field generation. Magnetic field is amplified only by compression. Then 
the shocked magnetic energy per particle is ^ few eV, because the shock has a compression factor 
~few. But the proton temperature in the shocked plasma is ~ few GeV, giving ^b ~ 10~^, which 
is at least 6 orders of magnitude smaller than required by the synchrotron phenomenological model 
of GRB afterglow. 



In this work, we assume that the synchrotron GRB model is correct. Then the magnetic field 
should be somehow generated by the shock. For all practical purposes, the unshocked plasma 
is unmagnetized. The nature of this shock-generated magnetic field is an interesting unsolved 
problem. Here we try to guess an answer, but we have to acknowledge the speculative nature 
of our discussion. The right way to solve this problem is direct numerical simulations of a three 
dimensional collisionless (Vlasov) plasma. 



2. Shock Dynamo 

That collisionless shocks can generate strong, that is near equipartition, magnetic fields is well 
known (eg. Sagdeev, 1966). Low-dimensional numerical simulations illustrating this process are 
also available (eg. Kazimura et al 1998, Appendix). The mechanism of the magnetic field generation 
is the Weibel instability, which can be explained as follows. Within the shock transition region, 
the distribution function of charged particles is strongly anisotropic. Particles move predominantly 
perpendicular to the front. These particles might be thought of as straight currents. Currents of 
the same sign attract each other, leading to bunching of electrons moving in one direction and 
protons moving in the opposite direction. This current bunching generates magnetic fields. 

The length and time scales characterizing the instability can be obtained from dimensional 
analysis. The instability is collisionless — it is described by the Vlasov system of equations (that 
is collisionless Boltzmann plus Maxwell). The Vlasov system contains the following dimensional 
quantities, nm, ne, and c. Here n is the unshocked plasma density, m is the proton mass, c is the 
speed of light.Q We will assume that the shock is mildly relativistic, so it is not characterized by 
any large or small dimensionless numbers. Then there is only one time scale in the problem, the 
inverse plasma frequency, ujp^, where ujp = (Airne'^/m)^''^. The only length scale is the skin depth, 
6 = c/ojp. 

Now the initial stage of the Weibel instability might be described as follows. The shock width 
is ~ 6. Within the shock, a near equipartition magnetic field is generated. The field isotropizes the 
distribution functions of charged particles. The characteristic length scale of the generated field is 
^^ 6. 

But what happens many deltas downstream (Gruzinov and Waxman, 1998)? The standard 
plasma-physical answer would be: (i) Strong electromagnetic fields occupy a layer several 5 wide. 
Within this layer, electromagnetic fields play a role of collisionality, bringing the plasma up the 
shock adiabatic. Many deltas downstream, there is only shocked plasma and virtually no fields. 
Electromagnetic fields are week, because they decay by phase space mixing (Landau damping) . We 



^For the sake of simplicity, we will exclude a small dimensionless number me/m, where me is the mass of electron, 
from the list of parameters. In reality this number might be important, for example, Sagdeev (1966) argues that 
non-relativistic Weibel saturates at ^b ~ me/m. 



have confirmed this magnetic field decay scenario, which is the common colhsionless shock wisdom, 
by two-dimensional numerical simulations (Appendix). 

But: the decaying field scenario is ruled out by astronomical observations, so long as we 
believe that GRBs and GRB afterglows are synchrotron emitting shocks. Synchrotron emission of 
a few-skin-deep layer would be negligible. 

The two most natural possibilities consistent with the synchrotron GRB model are: (ii) strong 
small-scale fields (scale ~ 6) exist downstream, (iii) strong large-scale fields (scale ~ /, /is proper 
distance from the shock) exist downstream. Now (ii) is, in fact, impossible. Magnetic fields don't 
live on length scales 6. They can be born there, and then they decay there on a time scale ~ uj~^ 
due to phase space mixing. Our numerical simulation confirms that magnetic fields really die out 
on skin depth scales (Appendix), at least in 2D. We will assume that the decay of microscopic, 
skin-depth magnetic fields by the phase space mixing is a true phenomenon, not an artifact of two 
dimensions. In 3D the phase space mixing should be even easier. 

Then we are left with the possibility (iii) - strong (but maybe ^b ^ w-e/^^-p) large-scale (length 
scale of order distance from the shock) magnetic fields are somehow generated by a colhsionless 
shock propagating into an initially unmagnetized plasma. We do not understand how this happens. 
We suspect that particle acceleration is an essential part of the process.^ The presence of an 
energetically important population of accelerated particles invalidates the MHD approximation. 
And one does have to go beyond MHD, because MHD shocks are stable and therefore do not 
generate magnetic fields. 

We think that the problem of magnetic field generation in shocks will be solved by direct 
numerical simulations. Our amateur numerical program seems to be sufficient for a 2D case. A 
carefully designed code should solve the 3D problem. The advocated 3D simulation of a colhsionless 
shock propagating into an unmagnetized medium should be performed regardless of our arguments.^ 
This simulation might prove our scenario wrong, but any result should be of great interest. Say, one 
discovers that magnetic fields are not produced. Then, what is the nature of the GRB afterglows? 



3. The First Magnetic Fields 

We have suggested that relativistic colhsionless shocks generate strong large scale magnetic 
fields. However we do not see how the relativistic effects might be important. We therefore 
hypothesize that non-relativistic shocks also generate strong large scale magnetic fields. 

This means, in particular, that the first cosmic magnetic fields could appear together with the 
first virialized objects. The magnetic fields are generated in virialization shocks, and they are born 

^In 2D, we see no magnetic field generation on large scales and no particle acceleration. 

''Such simulation is needed to understand particle acceleration. Postulating "scattering clouds" is not good enough. 



strong and large-scale. The primordial magnetic field theories are not needed to explain "the seed 
fields" for galactic dynamo. 



A. Weibel instability in 2D3V 

Here we show that magnetic field generated by a 2D3V (two-dimensional in space and three- 
dimensional in velocity space) Weibel instability decays due to phase space mixing after a few 
plasma times. We perform a simplified numerical simulation to illustrate this effect. 

We use dimensionless units e = m = c= 1. Electromagnetic field is 

E = {0,0, -dtA), B = {dyA,-d^A,0), (Al) 

where the z-component of the vector potential is ^ = A{x,y,t). The absence of the electrostatic 
field is explained below. Non-relativistic charges move in a 2D space, r = {x,y), with 3D velocities 
(v, m) = {vx,Vy,u). This means that charges are rigid rods elongated along the z-coordinate. The 
equations of motion for positive charges are 

r = V, V = uVA, u = -dtA - v • VA. (A2) 

To simplify our basic system, we will assume that at t=0, and therefore at all later times, each 
positive charge {r,v,u) is matched by a negative charge (r, v, —u). Under this symmetry, there is 
no electrostatic component of the field. The only field equation is 

- V^A + d'^A = j, j = n'^ f (fvdu fu. (A3) 

Here il^ = Svrne^/m, n is the density (the unperturbed distribution function / = ^(v, u) is 
normalized J d^vdu F = 1). 



A.l. Linear Theory 

To analyze the linear Weibel instability, we follow the standard Vlasov route (eg. Lifshitz & 
Pitaevskii 1981). We represent the equations of motion by the Boltzmann equation: 

5i/ + V • drf + uVA ■ d^f - {dtA + V • VA)duf = 0. (A4) 

For a perturbation oc e~*'^*+*'^'", the perturbed part of the distribution function is 

6f = A{{uJ-k■^r)-\k■^^ + ^u)F. (A5) 

We use this in (A3), and after integrating the current density by parts, we get the dispersion 
law (that is relation between to and k): 



where < ... > is the average over the unperturbed distribution function. 

If, as in our numerical simulation, at the initial time all velocities are parallel to z, (A6) gives 
the instability growth rate 

This means that the fastest growing modes are small-scale, with wavelength A ^ $7^^. The charac- 
teristic growth rate is ~ uQ. 

A. 2. Nonlinear Theory 

As the instability develops, magnetic field turns the particles' trajectories, and a non-zero 
velocity dispersion in the (x,y) plane appears. This stabilizes the short wavelengths first, because 
in the first non-vanishing order in v, (A6) gives, for k ^ Q, 

-co^^n^<u^> - ^ ^^Y^ fc^ (A8) 

2 < u2 > 

One expects maximal convertion of kinetic into magnetic energy at about the time when one 
Larmor circle in the generated magnetic field has been completed by a typical particle. This should 
occur at length scales somewhat larger than l/fi, because the generated field has to be somewhat 
subequipartition, that is less than Uu. 

Further evolution is governed by conservation of the mean squared potential. As follows from 
translation invariance along z, the generalized momentum, p = u{t) + A{x{t),y{t)) is conserved for 
each particle. In particular, / = ^p'^/N, where the sum is over all N particles, is conserved. / has 
dimensions of the squared potential, that is B^L^, where B is the magnetic field and L is the length. 
One may then suggest the following scenario for the magnetic field evolution. The length scale L 
grows, / is conserved, therefore the magnetic energy < B^ > decreases as I /L?. To determine L 
as a function of the evolution time t, we assume that L/t is close to the Alfven velocity, which is 
oc B. It follows that < B'^ Xx t^^. 



A. 3. Numerical Simulations 

In our simulations, we integrated the particle trajectories (A2), and solved the electrodynamics 
equations dtA = —E, dfE = —V'^A — j on a grid. Space was a periodic square of size 1. We used 
the simplest straightforward numerical algorithm. With a 100x100 grid for an electromagnetic field, 
we needed N=2,000,000 particles to compensate for the discreetness effects. For Q = 100, with the 
time step of 0.0005, energy is conserved to 0.3% , and the mean squared potential (see §A.2) is 
conserved to 3%. 



Initial conditions are as follows. Electromagnetic field is zero. Particles are at random spatial 
positions, with v = 0, and u randomly chosen from an interval (—0.2,0.2). This velocity might 
seem to be too high for a non-relativistic approximation to apply, since in our units the speed of 
light is 1. However, as we have checked, even the fastest particle in our run does not get faster 
than about 0.5. Results are shown in figures 1-3. 

At i = 1, 16% of energy is converted into magnetic energy. The mean magnetic field is -B ~ 3. 
Since velocity \sv ^ 0.1, the corresponding Larmor radius is ~ 0.03, close to the length scale seen in 
Fig.l. As seen from Fig. 2, the length scale of the field grows. As we have said, this is accompanied 
by the magnetic energy decay. As seen from Fig. 3, magnetic energy is indeed described by the t~^ 
law. 
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Fig. 1. — Generation of the magnetic field. Particles with positive u- velocities and the magnetic 
field at t = 1, when the field energy is 16% of the total energy. 




Fig. 2. — Decay of the magnetic field. Particles with positive u-velocities and the magnetic field at 
t = 10, when the field energy is 4% of the total energy. 
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Fig. 3. — Growth and decay of magnetic energy. Circles: numerical data. Solid line: a two- 
parameter fitted i~^ law: 0.57(2.7 + t)~^ 



